# Load necessary libraries and functions
# load functions
se <- function(x) sqrt(var(x, na.rm=T)/sum(!is.na(x)))

# load libraries
library(tibble)
library(lme4)
library(car)
library(ggplot2)
library(tidyr)
library(reshape2) 
library(patchwork)
library(knitr)
library(kableExtra)
library(maps)
library(dplyr)
library(sf)
library(rgdal)
library(officer)
library(rstatix)
library(gt)

